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Cellular cargo can be bound to cytoskeletal filaments by one or multiple active or passive molecular 
motors. Recent experiments have shown that the presence of auxiliary, nondriving motors, results 
in an enhanced processivity of the cargo, compared to the case of a single active motor alone. 
We model the observed cooperative transport process using a stochastic model that describes the 
dynamics of two molecular motors, an active one that moves cargo unidirectionally along a filament 
track and a passive one that acts as a tether. Analytical expressions obtained from our analysis 
are fit to experimental data to estimate the microscopic kinetic parameters of our model. Our 
analysis reveals two qualitatively distinct processivity-enhancing mechanisms: the passive tether 
can decrease the typical detachment rate of the active motor from the filament track or it can 
increase the corresponding reattachment rate. Our estimates unambiguously show that in the case 
of microtubular transport, a higher average run length arises mainly from the ability of the passive 
motor to keep the cargo close to the filament, enhancing the reattachment rate of an active kinesin 
motor that has recently detached. Instead, for myosin-driven transport along actin, the passive 
motor tightly tethers the cargo to the filament, suppressing the detachment rate of the active 
myosin. 



PACS numbers: 87.10.Mn, 87.15.hj, 87.16.Nn 

I. INTRODUCTION 

To carry out its functions, a living cell requires 
the precise spatiotemporal organization of many macro- 
molecules. Trafficking of molecules within the cytoplasm 
can be mediated by distinct processes including diffu- 
sion, polymerization, and active transport^. A variety of 
transport mechanisms may arise from the physical prop- 
erties of the diverse cargoes being transported. For ex- 
ample, in the case of large cargoes, such as organelles, 
mRNA or virus particles, diffusion may not be suffi- 
ciently fast nor be spatially controlled. These cargoes 
are often transported by motor proteins that processively 
move to and from the nucleus along specific cytoskeletal 
filaments^. 

The cytoskeleton is typically composed of three types 
of filaments: microfilaments (e.g. actin), microtubules 
(e.g. tubulin a and /?), and intermediate filaments (e.g. 
lamins) 3 . Molecular motors most often associate with 
and process cargo along actin and microtubules 4 . These 
two filament types are structurally very different from 
each other. Microtubules (MT) are thicker (25nm diam- 
eter) and have a specific radial orientation with respect to 
the cell nucleus. Actin filaments are more randomly dis- 
tributed near the peripher y of the cell, and are less thick 
(8nm diameter) than MTs^l Moreover, filaments are 
directional. The ends of a microtubule are structurally 
different and labelled "positive" or "negative," while the 
ends of actin filaments are "pointed" or "barbed." Ac- 
cordingly, there arc different types of motor proteins as- 
sociated with not only different filament types, but with 
the direction of transport along these filaments. The MT- 
specific motor proteins are kinesins (e.g., kinesin I, II) 



that can move along MTs away from the nucleus on the 
positive direction, and dyneins that move in the oppo- 
site direction towards the negative end of a microtubule^. 
Various forms of the actin-specific motor myosin trans- 
port cargo towar d th e barbed (e.g. myoV) or pointed 
(e.g. myoVI) endJ^. Since each motor is highly selec- 
tive, and cargoes need to be moved on both directions of 
each filament, there are other proteins/cofactors that fa- 
cilitate molecular transport by associating with specific 
motors and filaments. For example, dynactin is a co- 
factor of dynein that enhances both theprocessivity of 
dynein and its affinity to certain cargoes^. 

Single molecule imaging methods have been pivotal in 
the experimental study of molecular motor dynamics^'. 
Such advanced techniques have allowed researchers to 
dissect various aspects of molecular motor mediated 
transpori) 10 | 11 | 12 | 13 | 14 | 15 [ The identification of motor pro- 
teins, their structure a nd propert ies have also led to sev- 
eral theoretical studies ^ * 17 * 18 * 19 ! that have further im- 
proved our understanding of how a single motor protein 
is able to move a cargo along a cytoskeletal filament. 

Since many associated proteins/cofactors affect trans- 
port dynamics^, experimental and theoretical investiga- 
tions of model systems that include only one motor and 
the tracks on which they bind do not yield a complete 
description of molecular motor based transport in vivo. 
Therefore, other recent studies have focused on how co- 
operativity among different molecular motors can facili- 
tate cargo transport along st raight, branched, and inter- 
secting cyt oskele tal filament d 7 * 20 * 21 !. In the experiments 
of Ali et ai r 2 ' 23 ' the cargoes were fluoresccntly labeled 
quantum dots (Qdots)^ simultaneously attached to one 
kinesin and one myoV motor. While analyzing the dy- 
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FIG. 1: Schematic of movement along a filament track of a 
cargo (blue circle) with both an active (green) and passive 
(red) motors attached. Each figure represents one of three 
possible states the cargo complex can be in before detach- 
ment, together with the rate constants that lead the system 
out of a state. Note that the complex moves processively only 
if the active motor is attached. Only in state (0, 1) can the 
system diffuse. 

namics of myoV transport in the presence of both actin 
and MT filaments, the authors discovered that myoV, 
besides processing along actin, is also able to associate 
with, and diffuse along microtubules. In further work 23 , 
the same authors showed that when both myoV and ki- 
nesin motors are attached to a single cargo, the proces- 
sivity of the entire assembly is increased on both MT and 
actin filaments. 

In this paper, we develop a stochastic model for 
the cooperative enhancement of kinesin and myoV 
as discovered irPU. Related stochastic models, both 
discrete and continuous hav e bee n effectively used to 
model single motor dynamic d 18 * 25 1 and cooperativity^. 
When applied to single motors, these stochastic models 
describe the stepping dynamics of typically two-headed 
molecular motors (e.g. myoV or kinesin) that walk 
hand-over-hand along the filament, where the trailing 
and the leading heads exchange roles at each step. In 
our analysis, we take a coarser approach by treating the 
entire kinesin-myoV-cargo complex as a single structure 
that can exist in four possible states corresponding to 
different association states of each motor with each 
filament type. From the model we extract and briefly 
discuss the expressions for the mean run length and the 
first passage time before detachment. In the Analysis & 
Discussion section, we fit the model to the experimental 



dP a {x,t) | v dP a (x,t) = D d 2 P a (x,t) 
dt a dx a dx 2 



Here, v a is the velocity of the cargo in state a. This veloc- 



results fronPHlSl The fitting highlights a qualitative 
difference between cargo transport on actin and mi- 
crotubules. We also discuss the dependence upon the 
initial conditions of the system and perform a sensitivity 
analysis on the unknown parameters. 

II. STOCHASTIC MULTISTATE TRANSPORT 
MODEL 



Consider the molecular transport system of a cargo 
with two motor-proteins attached, one that acts as an 
active motor (e.g. kinesin on MT) and one that act as a 
passive motor or tether (e.g. myoV on MT). We denote 
the state of the engine a a and of the passive tether a p 
by an ordered pair a — (<J a ,a p ), where a aiP = 1 if the 
active/passive motor is attached to the filament track 
and where <j a ,p = otherwise. Using this notation, state 
a = (1,1) corresponds to the case where both engine 
and tether attach cargo to the filament. State (1, 0) rep- 
resents a cargo complex whose engine is attached, but 
where the tether has detached from the filament (al- 
though it remains attached to the cargo). Conversely, 
(0, 1) denotes the case where the engine has detached 
from the filament, but the tether still holds the cargo on 
the track. Finally, when both motors have detached from 
the filament, the complex reaches the (0, 0) state. The 
states of the cargo system are schematically shown in Fig. 
[I] We only consider states in which, at any given time, 
at most only one active motor and one passive motor 
tether can attach a cargo to a filament. This assumption 
is reasonable because the procedure used irP^ for motor 
attachment predicts that 95% of the Qdots used as cargo 
have only one active motor attached. While there are no 
predictions about the number of tethers attached to the 
cargo, the relatively small size of the Qdots used in the 
experiments (~ 15 nm diameter) suggests that simulta- 
neous attachment of an active motor and multiple tethers 
is highly unlikely. 

The cargo complex can move processively along a fila- 
ment only if the active motor is attached to both cargo 
and filament. In our notation, this corresponds to states 
with first index one, e.g. (1,1) or (1,0). In state (0,1), 
the cargo is either diffusing or immobile, depending on 
the property of the passive motor. Within this context, 
we can write the Master equation for the probability den- 
sity function P a (x, t) that the motor-cargo complex is in 
state a between position x and x + dx at time t: 



£ {W(a,a')P a ,(x,t) - W(a',a)P a (x,t)} . (1) 

a' 

i 

ity will depend on the specific properties of the motor- 
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TABLE I: Description of transition rates W(a, a') in Eq. [l 



State Transition 


Rate 


Description 


(1.1) 


- (0,1) 


d a 


Detachment rate of motor from 


motor-tether complex. 


(1,1) 


-(1,0) 


dp 


Detachment rate of tether from 


motor-tether complex. 


(0,1) 


— (1,1) 




Attachment rate of motor to tether 
only complex. 


(1,0) 


-(1,1) 


kp 


Attachment rate of tether to motor 
only complex. 


(1,0) 


-(0,0) 


Ma 


Detachment rate of motor from 


motor only complex. 


(0,1) 


-(0,0) 


Hp 


Detachment rate of tether from 


tether only complex. 



protein and on any externally applied forces. For ex- 
ample, as has be en we ll established under many experi- 
mental condition d 14 ! 21 ! opposing forces applied to motor- 
driven cargoes linearly decrease their velocity. 

Consistent with observations, we set the diffusion con- 
stant D a = in Eq. [TJwhcn an active, driving motor is 
attached, suppressing random diffusional motion. Con- 
versely, when only a passive tether is attached, the mo- 
tion of the cargo is Brownian and D a > 0. 

The last term in Eq. [I] represents transitions among 



sP(q,s)- 



where P(q, s) = P(x, s)e~ lqx dx. In the above equa- 
tion, and in the remainder of the paper, we use the sub- 
scripts a, p and ap to indicate the quantities and/or ex- 
pressions that are characteristic of a transport complex 
consisting of an active motor only (a), a passive motor 
only (p) or both (ap). In Eq. [5J we assumed that motion 
is purely convective when an active motor attaches the 
cargo to the filament, and therefore set D a = D ap = 0. 
As indicated by the experimental data irP^, the passive 
tether does not noticeably affect the transport velocity of 
an active motor. Therefore, we also set v ap — v ai where 
v ap is the intrinsic velocity of cargo in the state (1,1). 
Finally, since the passive motor acting as a simple tether 
cannot induce drift along a filament track, we set v p = 0. 
After imposing these physical constraints, we solve the 



binding states ex. The corresponding rates W(a, a') are 
assumed to be constant and are defined in Table |U Note 
that we make the physically reasonable assumption that 
a cargo complex in state (1,1) cannot have both motors 
detach simultaneously from the filament track, allowing 
no transition between state (1, 1) and (0, 0). 

We will analyze the model given by Eq. I 
by defining the probability density vector P(x,t) — 
(f(o,o), ^(0,1); -^(1,0), P(i,i)) T ■ If we assume that the fil- 
ament track is infinitely long and that the cargo is at 
position x — at initial time t — 0, the initial condition 
is P(x, t = 0) = (0, a, /3, 1 — a — (3) T S(x), where a is 
the probability that the cargo complex is initially bound 
to the filament only by the passive tether, and (3 is the 
probability that the cargo complex is initially bound only 
to the active motor. While it is experimentally difficult 
to quantify a and (3, it is relatively straightforward to 
determine how our main results for residence times and 
run lengths depend on the initial conditions. We can 
therefore establish, a posteriori, the importance of a, (3 
in the measured results. Thus, by studying how certain 
estimated quantities depend on the initial conditions, we 
can determine the significance and usefulness of experi- 
mentally pinpointing the exact values of a and (3. 

The analysis is facilitated by defining the Laplace 
transform in time P(x, s) — J °° P(x, t)e~ st dt, and by 
taking the dual Laplace-Fourier transform of Eq. [lj 



P, (2) 



I 

Master equation to obtain analytical expressions for the 
mean run length and time. We will use these expressions 
to compare our model with experimental data frorrPH. 



A. Mean Run Lengths 

Having more than one motor attached to a cargo com- 
plex typically results in an improved transport efficiency. 
Experiments by Ali et aZP^l show significant increases in 
the processivity of a cargo when it is also attached to a 
passive motor. Within our model, the measured proces- 
sivity is equivalent to the mean run length, (X ap ), of the 
cargo complex before detachment. To find its expression, 
we construct the detachment flux density: 



/ \ 



\l-a-p/ 



/0 \ 

D p q 2 

iqv a 

\ iqv a J 



( fi p /i a \ 

-fe a - fi p d a 

-fc p - dp 

V fe a fcp -rf a -dp J 
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TABLE II: Typical values of parameters and mathematical 
quantities for kinesin/myoV cargo transport. 





f\ /I non \/ q Inn 




rvei. 




1 0.76 fun 
1 1.7 /im 


myoV on actin 
kinesin on MT 


1231 


(X ap ) 


[ 1 DQ I ITT) 

J IiUl/ Lai lit 


mvnV-lf lnpQin on aptin 

111VU V lYllUjOlll KJ1.L GLV-'Ulll 


l23l 


I 3.7 fim 


myoV-kinesin on MT 




V a 


f 0.46 [im/s 


myoV on actin 


nu 


\ 0.88 fim/s 


kinesin on MT 




fJ-a 


| o.60 s _1 


tit^/ti \/ rrriim a p+~ i n 
iiiyu v ii win oii^Liii 


l23l 


\ 0.51 s" 1 


kinesin from MT 






> 0.2 s" 1 


kinesin to MT 


1231 




> 0.02 s" 1 


myoV from MT 


1351 


D 


0.11 - 0.26 fim 2 /s 


myoV on MT 


22 23 



Jap(x,t) = // a P(l, ) + Mp-P(0,l) 



(3) 



We can use the properties of the Laplace and Fourier 
transforms to find a compact expression for (X™), the 
771 th moment of the run length before detachment: 



J ap (x,t)dt 



Ma; 



dq 



Ma-P(l,0) 0, S = 0) + /JpP(o,l) S = 0) 

j ap (q,s = 0) 



9=0 



da; 



(4) 



By taking m = 1 we can find the expression for the mean 
run length of the motor-tether complex before detach- 
ment 



(X ap ) - 



(dp + fc p + /j a - /3/x a ) (fc a 4- /j p ) + d a /3/j p 



rfa f 1 p(kp 



Ma) 



<* [Mp ( fc P + Ma) + dpVp) 
dp A<a(fca 4- Hp) 



(5) 



As expected, Eq. |5]is independent of any diffusion of the 
passive tether since diffusion on average does not con- 
tribute to the mean displacement. The mean run length 
is a monotonically decreasing function of a for all the 
physically realistic (i.e. positive) values of the model pa- 
rameters. On the other hand, the dependence of (X ap ) 
on (3 will be monotonically increasing (decreasing) if the 
term <i a /i p — /i a (fc a 4- /J p ) is positive (negative). The pro- 
cessivity of a cargo complex initially in state (1,0) (when 
(3 = 1) is greater than that of a cargo beginning in state 
(1, 1) (when (3 — 0) only when the detachment rate d a 
of the active motor from the state when both motors are 
attached is greater than /J a (l + k a /fi p ), where /i a is the 
detachment rate of the active motor by itself. 

Establishing the dependency of Eq. [5]upon the remain- 
ing parameters does not lead to an easy expression unless 
we make some simplifications. We can focus on the limit 
in which all cargo complexes have both motors attached 
to the track (a = (3 = 0) at time zero. This is also the 
underlying assumption used irPH. In this case, we find 
that the mean run length before detachment reduces to 



(X a 



(hp 



Ma) (fca 4" Hp)v a 



Mprfa(fcp + /J a ) 4- /J a dp(fca 4- /Jp) ' 



(G) 



As expected, (X ap ) is a monotonically increasing function 
of k a , kp and v a , and a monotonically decreasing function 
of d a , dp, /J a and fi p . 



The probability density flux J a out of the state where 
the driving motor is attached can be expressed as 



J a (x,t) = /i a -P(l, )(M) + d a P(i,i)(x,t). 



(7) 



Since the final two states (0, a p ) are both absorbing, we 
must also set k a — in the evaluation of J a (x,t). The 
expression for the moments of this density flux is analo- 
gous to Eq. [4] The mean run length of a cargo complex 
conditioned on being bound by an active motor is given 

by 



(Xa) 



(3 (d a - /i a ) 4- (dp + k p + /x a ) 
rfpMa 4- d a (k p + /j a ) 



(8) 



and is a monotonically increasing (decreasing) function 
of /3 if d a > /i a (d a < /i a ). Similarly, a simple analysis 
of Eq.[8] reveals that if (3 — the mean run length (X a ) 
in state (l,a p ) is a monotonically decreasing function 
of fc p and a monotonically increasing function of d p if 
d a > Ma- If d a < /J a , the mean run length (X a ) in state 
(1, <t p ) is a monotonically increasing function of fc p and 
a monotonically decreasing function of d p . When the 
detachment rate of the active motor is independent of its 
state (e.g., /J a = d a ), (X a ) = v a /fi a is independent of the 
initial conditions and the properties of the tether. 

Finally, the mean run length of a cargo bound only by 
a passive tether (state (0, 1)) can be found from 
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(0,0) 




time 



FIG. 2: A representative trajectory. The processive runs dur- 
ing states (1, 1), (0, 1) are assumed to have the same velocity 
independently of the state they are in. The depicted cargo 
movement is characterized by three processive runs, all with 
the same slope (e.g. same velocity) and three diffusive events 
due to the cargo complex being in state (0, 1) before ccomplete 
detachment at X. The dashed red line represents a straight 
trajectory that defines an effective velocity of a complete run 
before detachment. 



J p (x,t) = (/i p + fc a )P(o,i)(a;,*), 



(9) 



and setting the rates into (0,1) from all absorbing states 
d a = (3 — 0, yielding (X p ) = 0. This trivial result stems 
from the drift-free nature (v p = 0) of the (0, 1) state. 



B. First Passage Times 

Together with the mean run length, the data irPS 
include the average velocities for all experimental mo- 
tor/tether configurations. To use this data for parameter 
fitting, we derive analytical expressions for the average 
duration of a run, i.e., the first passage time to the de- 
tached state. In order to evaluate the detachment time 
of the driving motor anywhere along the filament track, 



we consider the conditional mean first passage time T a 
out of states (1, a p ), for which a = 0, as follows 



io°° di -Go da; j «0m) 

1 dJ a (q = 0,s) 



J a (q = 0, s) 



ds 



(10) 



Upon using Eq.[l0] and fc a = together with the initial 
condition P(x, 0) = (0, 0, /?, 1 — (3) T 5(x) we find 



T 



P(d a - /l a ) + (d p + fcp + /i a ) 



We can use this result in conjunction with the expression 
for the mean run length from Eq. [8]to obtain an estimate 
of the velocity V a of the cargo conditioned on it being 
attached by the active motor: 



(Xa) 



(12) 



Since the passive motor does not affect the motion of the 
cargo, the velocity when an active motor is attached is 
independent of the binding and unbinding of the passive 
tether. This velocity is that within a single processive 
run. The velocity averaged over a trajectory composed 
of both processive and diffusive runs will typically be 
smaller, as discussed below. 

We now compute T ap , the mean first time to detach- 
ment (state (0,0)). Using analogous notation, we find 



T, 



ap 



J^dtf^dx t.J ap (x,t) 

1 dJ ap (q = 0,s) 



J ap (q = 0,,s) 
from which we obtain 



ds 



(13) 



s=0 



T = 

1 ap 



(4 + fc a ) (fcp + A^a ~ PfJ-a) + (/3«q ~ 

da, (fcp 



jU a ) + (1 - oQ(fc p + aQVp + d p (k a + a(^ a - M P ) + M P ) 

+ A* a )/X p + (ip^a(fca + Mp) 



(14) 



As in section II A we can study the dependence of (T ap ) 
on the kinetic parameters by considering the a = f3 = 
limit: 



d v (K + ^ p ) + (fcp + ^ a )(*4 + fca + Mp) 
da(fcp + A*a)Aip + dp^ a (fc a + A*p) 



(15) 



This expression is monotonically decreasing with respect 
to /z a an d Mpi an d can be used in conjunction with Eq. 
[6] to obtain an expression for the mean velocity of cargo 
transport in presence of an active and a passive motor 
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A. Parameter Fitting 








0.25 0.5 



0.75 1 



1.25 



FIG. 3: Surface plot of the parameter space that satisfies the 
conditions of Sec. III.A.l for the mean run length and first 
passage time. While both k p and d p can have values along 
the positive real line, fi p has an upperbound. If /j, p > 1.25s~ 
there are no physically realistic values of k p and d p that satisfy 
the conditions {X ap ) = 3.7 and T ap = 5.0. For this plot, 
v a = 0.88 (um, jU a = d a = 0.52 s~\ k a = 1.48 s -1 . 



(Va P ) = 



d p (k a + n P ) + (kp + ^ a )(cZ a + fc a + n P ) 



< v a . 

This result indicates that although in states (1,1) and 
(1,0) the cargo drifts with velocity v a , the entire run con- 
sists of alternating phases of drifting and diffusive states, 
resulting in an average effective velocity V ap < v a , with 
the two velocities being equal only under certain regimes 
such as for slow active motor dissociation from state (1,1) 
(e.g. 4 < 1) or for very fast tether dissociation from 
state (0,1) (e.g. fi p 3> 1). From the above expression 
one can estimate the ratio x °f convection times to dif- 
fusion times as follows 



\V ap -V a \ 



d a (k p + fx a ) 



dp(fc a + Hp) + ( k P + Ma)(<4 + fc a + (JLp) 



The above algebraic expressions for the mean length 
(e.g. Eq. [5]) and duration (e.g. Eq. 14 1 of a run out 
of the various states constitute our main mathematical 
results. 



III. ANALYSIS & DISCUSSION 

Having obtained analytical expressions for the mean 
run length and duration, we can use experimental data 
frorrP^l to study parameter assignment within our model. 
Because the experimental data refer cooperative motor 
transport along microtubules and actin, we will also di- 
vide the analysis into the two corresponding cases. 



In the absence of external forces v a is constant, 
v p = 0, and for fixed initial conditions, the model is 
characterized by 8 parameters. Only for a few of them 
it is possible to extract estimates from the available 
literature (see Table |n|. However, we can use certain 
biophysical constraints stemming fronP2l to reduce 
the parameter space as much as possible. Throughout 
this subsection the analysis is performed with a = (3 = 0. 

III.A.l myoV-kinesin transport on MT. Experi- 
mental results fronP^I show that the presence of a passive 
myoV motor, in addition to an active kinesin motor in- 
creases the typical run length of the cargo by two-fold and 
slightly decreases the velocity by ~ 15%. The same data 
show that the velocity and mean run length of cargo in 
state (1,0) and state (1, 1) are essentially the same. Fig. 
[2] is a graphical representation of a possible cargo trajec- 
tory showing three processive and three diffusive runs. 
Since all processive runs are observed to occur with the 
same velocity, the presence of the passive myosin does not 
affect the drive of the active kinesin motor. Therefore, 
both states (1,0) and (1,1) are indistinguishable within 
each processive run. These observations suggest that we 
can assume d a = Ha in Eq. [I] We used this assump- 
tion and the values of v a , /i a from Table Hi to solve the 
system consisting of Eq. [6] and Eq. [15 with the con- 
straints (X ap ) — 3.7 fim and T ap = 5.0 s obtained from 
the experimental results in 23 for microtubular transport. 

The solution of this system leads to a specific value 
of k a — 1.48 s _1 , implying an average diffusion time of 
1/fca = 0.68 s, consistent with the experimental results 
irP21 l n fact, the average run length of a kinesin-myoV 
cargo complex on microtubules lasts about 5 seconds and 
covers twice the distance of a Qdot with only a single ki- 
nesin motor attached to it. Therefore, the typical cargo 
movement due to kinesin/myoV motor consists of two 
processive steps (needing ~ 2T a — 4.4 s, see^) and one 
or two diffusive ones before detachment since we assumed 
that the cargo is initially in state (1,1). Given such ob- 
servations, our obtained value of fc a is consistent with the 
experimental data. More specifically, it suggests that on 
average there will be only one single diffusive event in 
between two processive runs. 

For the three remaining parameters we find that as 
long as /ip < 1.25s , we can always find k p and d p 
that satisfy the physical constraints (X ap ) = 3.7 /im and 
T ap = 5.0 s. These results are shown in Fig. [3] The 
upper limit for /i p implies an average diffusion time of 
about 0.8 seconds while in state (0, 1) and before detach- 
ment. However, the experimental results iri 22 l 23 l show an 
average diffusion time between 40 and 60 seconds in the 
absence of kinesin. The observed diffusion times are con- 
sistent with the value for k a obtained above and suggest 
that the detachment of the cargo from the microtubule is 
most likely to happen while the motor is in state (1,0). 

Overall, these results suggest that myoV increases the 



FIG. 4: Surface plot of the values of k p , d p , and d a that satisfy 
the mean run length expression in Eq. |18| when, v a = 0.46 [im 
and /x a = 0.60 s -1 . For k p <C d p , d a tends to the value 
0.42s - . This value is equivalent to having a cargo complex 
consisting of only one active and no passive motor. 



processivity of the cargo complex by keeping kinesin 
close enough to the track so that its reattach is acceler- 
ated. The tethering by myoV occurs without reduction 
in the intrinsic velocity of kinesin. 

III. A. 2 myoV-kinesin transport on actin. The most 
striking finding observed using this experimental system 
is that the processivity is increased by the presence of a 
passive kinesin tether, but the average velocity remains 
unchanged, suggesting that the passive tether helps keep 
the active myoV attached to the filament track with- 
out affecting its velocity. In these experiments, one ob- 
serves longer, uninterrupted processive transport, with 
rare punctuated moments of diffusive cargo motion. This 
suggests that the system is predominantly in states (1, j) 
and that once the active motor detaches, the whole cargo 
system does too. These observations are consistent with 
the structural/molecular attributes of this system, since 
the electrostatic forces between actin and kinesin are too 
weak to significantly reduce the myosin driven cargo ve- 
locity. Moreover, since actin filaments are thin, once 
myosin detaches, the detachment of kinesin is also fast. 
But if myosin holds kinesin proximal to the actin fila- 
ment, the kinesin attachment rate is also fast, since free 
diffusion is hindered by the tether. Within this context 
d& 7^ Ma- We can also assume that fc a = and that 
fXp — ► oo j leading to the following expressions for the 
mean run length and first passage time 



(X, 



(k p + d p + fx a )i 



d p fi a 



■{Ta; 



Ma 



d a {kp + /ia) + rfpMa 

(18) 

The above expressions are the same as Eq.[8] and Eq.[TT] 
for the particular choice of initial conditions that we have 
used throughout this section. Moreover, from Eq. [12] 
we know that the velocity of this system is constantly 
equal to v a . Therefore, the model is able to predict 
the observed experimental behavior. We can now use 



FIG. 5: Plot of the dependence of k p , d p , and (i p on /3 for 
transport along microtubules. The values of the other pa- 
rameters are the same as the ones in Fig. [3] Increasing the 
number of cargo complexes in state (1,0) at time t — leads 
to smaller values of tether detachment from state (1, 1) for all 
values of k p and /i p . 



(X a p) — 1.09 \im and T ap — 2.59 s (obtained from Table 
1 irP3) to plot the parameter space that satisfies either 



one of the expressions in Eq. 18 but not both since they 



are redundant under the assumption (V ap ) = v a . Using 
the expression for the mean run length in Eq. |18[ we 
obtain the plot shown in Fig. [4] Here, we see that for 
k p ^> dp, d a reaches a maximum value of d a — 0.42 s _1 . 
This quantity is smaller than /i a , a result that confirms 
our interpretation that the increase in processivity is due 
to the tethers ability to prevent detachment of the active 
motor. 

From the properties of actin transport discussed thus 
far, it seems natural to consider an "effective" detach- 
ment rate /j, e e that captures the dynamics of cargo trans- 
port in this case: 



Meflf = /(Ma, Mp, kpi da, dp, V a ) 



(19) 



From the results irPH, W e obtain Meff ~ 0-42 s _1 . This 
value is the same as the maximum value obtained above, 
implying that the overall effect of kinesin is to prevent 
myoV from detaching from the actin filament, without 
affecting the intrinsic myoV transport velocity. 



B. Dependence on Initial Conditions 

We now investigate how the model fits data as a 
and f3 in Eq. [2] vary between zero and one. Based on 
the experimental methodology oP^, we expect both a 
and (3 to be small. In fact, only Qdots that show a 
displacement greater than 0.3 /im where included in 
the data, with shorter trajectories discarded, essentially 
reducing a to zero. The argument for (3 being small also 
relies upon the experimental methods of Ali et al. To 
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ensure that each Qdot had at most one active motor 
and at least a passive one attached to it, they prepared 
a solution where passive motors were in excess, with 
molar ratio 16:1^. 

III.B.l myoV-kinesin transport on MT. As dis- 
cussed above, in the case of microtubular transport, we 
can use the approximation d a — /i a . Unlike the case 
with simple initial conditions (a = f3 = 0), including 
this constraint in Eqs.[5| and [14| does not yield a simple 
analytical solution for the mean run length and first 
passage time before detachment. Therefore, we per- 
formed a numerical study of the dependence on the 
initial conditions and found that for every choice of a 
and (3 there is only one value of k a that satisfies the 
experimental results (X ap ) — 3.7 /j,m and T ap — 5 s. In 
addition, this value depends only on the initial fraction 
of cargoes in state (0, 1) and not on (3. In particular, 
k a is a linear function of a, with slope m — 1.26 s , 
giving us a range of predicted values from k a — 1.48 s~ l 
(if a = 0) to k a = 2.74 s" 1 (if a = 1). The increase of 
predicted fc a with a is not surprising since the higher the 
probability the systems starts in state (0,1), the faster 
the kinesin motor will have to bind to the microtubule to 
satisfy the given time constraint. Conversely, the values 
of the other free parameters in the model (i.e. d p , /i p , 
and fc p ) depend only on the value of (3. This dependence 
is shown in Fig. [5] for the limit cases (3 = and (3 = 1 
(in both cases a = 0). From this plot we notice how 
an increase in the percentage of cargo complexes in 
state (0, 1) at t = 0, shifts the parameter surface down 
along the d p axis. As a result the parameter space 
itself in Fig. [5] is reduced, since any combination of 
parameters resulting in rf p < is unphysical. From this 
analysis it seems that all possible initial conditions can 
explain the experimental data. None of the qualitative 
observations made in Sec. III.A.l would change, unless 
a > 0. In this case, the average cargo movement would 
consist of two diffusive steps, one at the beginning of 
the motion and one at the end of the first processive step. 

III.B.2 myoV-kinesin transport on actin. Using 
the same assumptions discussed in Section III. A. 2, we 
find the following simplified expression, 

, y v _ (f3(d a - + (1 - a)(k p + fe)) + d p (l - a) 

a a (Kp + Ma) + «p + /"a 

(20) 

indicating a mean run length that decreases linearly as a 
increases. This is physically expected since we assumed 
k a = 0, which implies that diffusional states are not al- 
lowed to transition to the processive (1,1) state. How 
{X ap } above scales with (3 depends instead on the differ- 
ence between \x a and d a . The detachment rate of myoV 
from actin in state (1,0) is about 0.6 s _1 , and we see 
from Fig. [4] that this value is never reached (we also ver- 
ified this result asymptotically). Then, the mean run 
length of the cargo complex is a decreasing function of 




FIG. 6: Plot of the dependence of k p , d p , and d a on /3 for 
transport along actin. The values of the other parameters are 
the same as the ones in Fig. [4] Increasing the number of cargo 
complexes in state (1, 0) at time t = leads to smaller values 
of active motor detachment rate from state (1, 1). Moreover, 
some of the values of k p and d p lead to negative values of d a 
although the limiting behavior for k p <C d p is still the same 
as for the case /3 = 0. 

the probability of being in state (1,0) at time t = 0. 

The overall effect of a higher value of (3 is to lower the 
best-fit value of d a and to reduce the range of physically 
meaningful parameters, as shown in Fig. [6] Similar to 
what we mentioned in Section III.B.l, the experimental 
data can be fit to all possible initial conditions. How- 
ever, if dp is too small, the condition for the mean run 
length frorrPSl cannot be satisfied if (3 is too large. This 
is reflected by the negative values of the (3 = 0.5 fit to 
e? a in Fig. [6] These results reinforce the notion that the 
increased processivity of myoV due to the presence of ki- 
nesin depends on the latter's ability to keep the active 
motor attached to the track for a longer period of time. 

C. Sensitivity Analysis 

We conclude this section by performing both local and 
global sensitivity analysis of our model output on the 
model parameters. Since our goal is to determine the 
effect of cooperation among different molecular motors 
on the processivity of cargo transport, we select (X ap ) 
(e.g. the mean run length before detachment) as the 
output of interest. 

The simplest local sensitivity analysis evaluates the 
partial derivatives of the mean run length before detach- 
ment with respect to each of the unknown parameters 
(e.g. d(X ap ) /ddp). This determines how sensitive the 
output is to quantitative changes in each of the kinetic 
parameters^. 

Since local sensitivity analysis is best suited to evalu- 
ating output that is linear in the parameters, we will also 
consider global sensitivity analysis on (X ap ). This analy- 
sis will determine which among the unknown parameters 
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A ft local sensitivity b * local sensitivity c d local sensitivity 




TABLE III: First and second order sensitivity indexes for mi- 
crotubule and actin cargo transport. To determine these in- 
dexes we sampled d p and k p uniformly in [0, 20] with step 
0.05. For the microtubule case we sample /x p in [0, 1.25] with 
step 0.01. For the actin case we sampled d a in [0,0.42] with 
step 0.01. 
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0.06 
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0.35 


0.82 
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0.1 


0.33 




0.23 




0.45 




0.37 0.81 



FIG. 7: Local sensitivity analysis of free model parame- 
ters in the case of kinesin/myoV cargo transport on micro- 
tubules. The parameters for all three of these graphs are 
v a = 0.88 fxm/s, d a = Ma — 0-52 s _1 , fc a = 1.48 s _1 , and 
d p = 0.5 s _1 . A: Sensitivity analysis with respect to fi p , 
in this case the z-axis represents d(X a p}/dfj, p . B: Sensitivity 
analysis with respect to k p , in this case the z-axis represents 
d(X ap ) /dk p C Sensitivity analysis with respect to d p , in this 
case the z-axis represents d(X ap ) / dd p . 



of the model would be most responsible for experimen- 
tal variation of the outputP3. This analysis is global in 
the sense that it spans all of the parameter space. It is 
model- free, and gives the same result as local sensitivity 
analysis if the analyzed models are linear^. Let us write 
^ = A*ai d pl d a , k p , k a , v a } as the input space, and f^, 
i = 1, 2, . . . , 7 as i-th input in Q. Then we can define the 
first-order sensitivity for a fixed input as: 



V[E{{X ap )\tt t )] 
V[(X ap )} 



(21) 



where E((X ap )\fli) is the expected value of the mean 
run length obtained by uniformly sampling over all other 
parameters Clj^ti, V[E((X ap )\Sli)] is the variance of the 
expected mean run length over the parameter f2i, and 
V[(X ap )] is the unconditional variance of the mean run 
length. The parameter with highest first order sensitivity 
index is the one which most influences the variation of 
the mean run length according to the global sensitivity 
analysis approach. Global sensitivity analysis can also 
be used to assess the joint effect of more than one input. 
We define the second order sensitivity (also known as 
two-way interaction) as 



V[(X ap )] 



(22) 



where E((X ap )\Qi, flj) is the expected value of the mean 
run length given fixed values of £li and £lj . Higher order 
global sensitivity indexes can be analogously denned. We 
apply local and first and second order global sensitivity 



analyses to both experimental cases: 

III.C.l myoV-kinesin transport on MT. Represen- 
tative results of the local sensitivity analysis for the mi- 
crotubule case are plotted in Fig. [7J Under certain 
regimes, /z p has the greatest influence on the mean run 
length before detachment, followed by d p , with k p as the 
least influential among the three parameters. To further 
investigate the results from the local sensitivity analysis, 
we determine the first and second order global sensitivity 
indexes for all free parameters. These results are listed 



in the first row of Table III We find that the parameter 



that is responsible for most of the variation in (X ap ) is 
d p . 

Both of the analyses predict that k p has the least 
influence on the mean run length, but they differ in 
their ranking of /i p and d p . This difference arises from 
the nonlinearity of Eq. [5] and the intrinsic differences 
among the two types analyses. Local sensitivity analysis 
suggests that if we could control the values of the 
parameters of the system, we would effect the largest 
changes in (X ap ) by altering the rate of detachment of 
myoV while in state (0,1). If experimentally, we are 
sampling parameter space, global sensitivity analysis 
predicts that by correctly determining d p we can achieve 
the most reduction in the variability of the mean run 
length. 

III.C.2 myoV-kinesin transport on actin. The lo- 
cal and global sensitivity analyses also give different re- 
sults in the case of cargo transport along actin filaments. 
Three representative plots of the local sensitivity analy- 
sis are shown in Fig. [8] These show that the detachment 
rate of myoV from actin when both motors are attached is 
the most influential parameter with respect to the mean 
run length before detachment. However, the first order 
sensitivity index of d a is about 1/3 smaller than the same 



index for d p (see Table III I, making the detachment rate 



of kinesin when in state (1,1) the parameter more re- 
sponsible for the variance in (X ap ). In this case, the 
sensitivity to d p is consistent with the role of the passive 
tether in preventing myoV detachment. 

Interestingly, both cases (actin and microtubule) have 
dp as the free parameter with highest first order sensi- 



10 



A d local sensitivity B k local sensitivity c d local sensitivity 




10 



k 1 \/ 0.5 

P d ; 



FIG. 8: Local sensitivity analysis of free model parameters 
in the case of myoV/kinesin cargo transport on actin. The 
parameters for all three of these graphs are v a = 0.46 /im/s, 



/ia = 0.6 s 



s , /j p < 1 and d p = 0.5 s 



Sensitivity analysis with respect to fj. p , in this case the z-axis 
represents d{X ap ) /<9d a - B: Sensitivity analysis with respect to 
fc p , in this case the z-axis represents d{X ap ) /dk p C Sensitivity 
analysis with respect to d p , in this case the z-axis represents 

d(x ap )/dd p . 



tivity and k p as the lowest one. Moreover, if we want to 
reduce by at least 80% the uncertainty in the determi- 
nation of the mean run length, we can do so by exactly 
determining d p and /i p for the mic rotu bule case and d p 
and d a for the actin case (cf. Table III). 



IV. SUMMARY & CONCLUSIONS 

We presented a stochastic model that describes the 
cooperative behavior between two different motors at- 
taching cargoes to a cytoskclctal filament. Of the two 
motors, one acts as engine, moving the cargo unidirec- 
tionally along the filament, while the other acts as a 
tether. We applied the model to the data frorrPSl, where 
the authors studied in vitro the cooperative behavior of 
kinesin and myosin V when moving cargo along micro- 
tubules and actin. Experimental visualization indicates 



significant diffusive dynamics for cargo transport along 
microtubules in the presence of myosin V, while cargoes 
transported along actin did not exhibit diffusive dynam- 
ics. A consistent interpretation of these observations is 
that along microtubules, tethers (myosin) predominantly 
enhances reattachment of the active motor (kinesin) . In 
the case of transport along actin, the kinesin tether acts 
to prevent detachment of the active myosin motor. This 
interpretation has been verified within our model. In 
fact, for the case of microtubule transport we found that 
the reattachment rate (fc a ) of kinesin to the filament when 
the tether is attached is three times faster than its cor- 
responding detachment rate. For the case of actin, we 
found that the detachment rate of myoV when kinesin is 
attached to the filament is slower than in the absence of 
the tether for all values in the explored parameter space. 

Although most of the modeling effort described in the 
literature use a very detailed representation of the physi- 
cal and chemical properties of the molecular motor under 
investigation, such analysis becomes much more complex 
when considering multiple cooperating motors. Coarser 
models such as ours can be straightforwardly extended 
to incorporate cargo systems with M active motors and 
N tethers. The state of such cargo complex will be char- 
acterized by a M + A-tuple, whose first M components 
have value 1 or depending on whether the active motors 
they represent are attached to the cytoskeletal filament or 
not. Conversely, the last N components will have value 1 
or depending on the attachment status of each passive 
motor. Proceeding as we did in Section [TTJ we can define 
a Master equation from which mean run lengths and first 
passage times to full detachment may be determined. As 
done in the present case, such a model could be used to 
find the optimal configuration of motor and tethers that 
fits the available experimental data. 
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